Singular reduction of implicit Hamiltonian systems 



Guido Blankenstein 
Department of Mechanical Engineering 
Katholieke Universiteit Leuven 
Celestijnenlaan 300 B 
B-3001 Leuven (Heverlee) 
£^ ■ Belgium 

o 

(N ; Tudor S. Ratiu 

H ■ Institut de Mathematiques Bernoulli 

CZ ■ 

Ecole Polytechnique Federate de Lausanne 
MA-Ecublens, CH-1015 Lausanne 
Switzerland 

\ Abstract 

Q 

This paper develops the reduction theory of implicit Hamiltonian systems admitting a sym- 
'"^ , metry group at a singular value of the momentum map. The results naturally extend those 

known for (explicit) Hamiltonian systems described by Poisson brackets. 
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1 Introduction 



Consider a symplectic manifold (M, oS) admitting a symmetry Lie group G (with Lie algebra denoted 
by 0) acting freely and properly on M, together with a corresponding equivariant momentum map 
P : M — > g*. In 19 it is shown that, at a regular value \i G Q* of the momentum map, the symplectic 
structure on M naturally reduces to a symplectic structure uj^ on the reduced manifold M M = 
P^ 1 (fi)/G fM , where G M is the coadjoint isotropy subgroup of G. Furthermore, the integral curves 
of a Hamiltonian vector field corresponding to a G-invariant Hamiltonian H € C°°(M) G project to 
integral curves of the reduced Hamiltonian vector field corresponding to the reduced Hamiltonian 
Hn € C°°{M^). This theory has been generalized in ^7] to the case of Poisson manifolds: at a 
regular value fi of the momentum map, the Poisson bracket {•, •} : C°°(M) x C°°(M) — > C°°(M) 
on M descends to a Poisson bracket {•>•}« on the reduced phase space ¥». Again, the Hamiltonian 
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flow defined by a G-invariant Hamiltonian H € C°°(M) reduces to a Hamiltonian flow on 
corresponding to £ C oc (M fi ). We refer to [Tl 1181 122] for some excellent overviews of this theory. 

Although by Sard's theorem regular reduction describes the "generic" case (i.e., the regular val- 
ues of the momentum map are dense in g*), interesting dynamics, such as bifurcation phenomena, 
occur at the singular values of the momentum map. Already in simple examples like the spherical 
pendulum moving with angular momentum zero (i.e., moving in a plane) and the Lagrange top 
(reducing the gravitational S 1 symmetry after regular reduction of the internal S 1 symmetry cor- 
responding to the homogeneous mass distribution), one recognizes the need for the investigation of 
these special cases. This task has been taken up in jHJ 01 EEH 1221 EH HZ] leading to the theory 
of so-called singular reduction of (symplectic and Poisson) Hamiltonian systems. See |§] for a nice 
overview and some worked out examples (including the spherical pendulum and the Lagrange top). 
The main difference with the regular reduction theory is that in the case of a singular value fi of the 
momentum map (i.e., the derivative, or tangent map, of P at points of P~ l (n) is not surjective), the 
reduced space is not a manifold. Therefore symplectic forms and (Hamiltonian) vector fields 
are not defined and, as a consequence, the reduced "Hamiltonian dynamics" cannot be written as a 
system of ordinary differential equations on (as in the regular case). However, since is still 
a topological space (relative to the natural quotient topology), a reduced Poisson bracket {-,-}^ 
on the space of (Whitney) smooth functions C°°{M lJL ) can still be defined. This bracket induces a 
Hamiltonian formalism that allows one to write the reduced Hamiltonian dynamics on the singular 
reduced space M^. The Hamiltonian flow corresponding to this Hamiltonian dynamics is exactly 
the projection of the regular Hamiltonian flow on M. Finally, in [1J 03 I1UI 111! 1241 127j it has been 
shown that the singular reduced space (resulting from a symplectic manifold (M, u)) may be 
stratified by symplectic manifolds, called pieces. The stratification is by orbit type decomposition. 
The Hamiltonian flow on leaves these pieces invariant and restricts to a (regular) Hamiltonian 
flow on each of the pieces. 

Recently, E] have generalized the regular reduction theory for explicit Hamiltonian systems 
to a regular reduction theory for implicit Hamiltonian systems, extending preliminary work in 
[HI EHj. Analogous to the symplectic form or the Poisson bracket in the classical theory, the 
underlying geometric structure of an implicit Hamiltonian system is that of a Dirac structure, 
defined as a maximally isotropic smooth vector subbundle of TM © T*M. This structure allows 
to define a Hamiltonian formalism generalizing the classical symplectic and Poisson formalisms by 
including the description of algebraic equations. That is, the "Hamiltonian dynamics" corresponding 
to a Dirac structure and a function H € C°°(M) consists of a set of differential and algebraic 
equations (and is therefore called implicit). Perhaps the most striking example of an implicit 
Hamiltonian system is that of a Hamiltonian system defined by a Poisson bracket on M, restricted 
to a submanifold of M that is not the level set of a Casimir function (this example actually motivated 
the definition of a Dirac structure in [S] ) . Using the notion of a Dirac structure introduced in [HI El , 
implicit Hamiltonian systems were defined in [SD E2 EH1 1201 El and successfully employed in the 
context of network modeling of energy conserving physical systems such as mechanical systems with 
(non)holonomic kinematic constraints, electrical LC circuits, electromechanical systems. We refer 
to [2J5] and references therein for more information; see also [5] for a detailed historical account. 
(Note that the Hamiltonian systems defined in are n °t truly implicit, but rather explicit 
systems.) The study of symmetries and reduction of implicit Hamiltonian systems evolved from 
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preliminary results in jSJ 1131 128j and led to a reduction theory for this class of systems described in 
[SJ|f)j. There it was shown that a Dirac structure on M, admitting a symmetry Lie group G with 
corresponding equivariant momentum map P, reduces to a Dirac structure on the reduced manifold 
Mn if fj, is a regular value of the momentum map. Furthermore, the (projectable) integral curves of 
the implicit system corresponding to a G-invariant function H £ C°°(M) G project to integral curves 
of the reduced implicit Hamiltonian system defined by the reduced Dirac structure and the reduced 
Hamiltonian H^. The theory generalizes the classical regular reduction theory for symplectic and 
Poisson Hamiltonian systems, as well as the recently developed reduction theories for constrained 
mechanical systems. Section briefly recalls the main results of EJ; for a discussion of the 
connection with the reduction theory for constrained mechanical systems we refer to El • 

The goal of this paper is to develop reduction theory for implicit Hamiltonian systems at singular 
values of the momentum map. We restrict our attention to a fairly general class of Dirac structures, 
described by a generalized Poisson bracket and a distribution of derivations (i.e., vector fields) on 
smooth functions. We consider the special subclass of symmetries that are symmetries of both the 
generalized Poisson bracket and the distribution. Using these ingredients, we prove that one can 
define a so-called topological Dirac structure on the singular reduced space (where for easiness 
of exposition we will take fj, = 0), representing the reduced Dirac structure. This topological 
Dirac structure, whose construction implicitly uses the Sikorski differential spaces (see |11| 26 ), 
defines a Hamiltonian formalism on the reduced space and the dynamics corresponding to 
an implicit Hamiltonian system with Hamiltonian H E C 00 (M fl ) are described. It is shown that 
the (projectable) integral curves of the implicit Hamiltonian system on M corresponding to a 
G- invariant function H € C°°(M) G project to "integral curves" of the Hamiltonian dynamics 
defined on the singular reduced space. The orbit type decomposition defines a stratification of the 
singular reduced space. It is shown that, under certain conditions, the topological Dirac structure 
restricts to regular reduced Dirac structures on the pieces (which are always manifolds) and that 
the Hamiltonian flow restricts to a regular Hamiltonian flow on each of these pieces. 

The paper is organized as follows. Section |2] gives a brief introduction to Dirac structures and 
implicit Hamiltonian systems. The basic results concerning symmetries and regular reduction of 
implicit Hamiltonian systems are recalled in section El Section 0] describes in a purely topological 
way how to reduce an implicit Hamiltonian system admitting a symmetry group. The result is a 
so-called topological Dirac structure on the singular reduced space (in general, not a manifold). It 
is shown that if the symmetry group acts regularly and the value of the momentum map is regular, 
and hence the reduced space is a smooth manifold, then the singular reduced implicit Hamiltonian 
system equals the regular reduced implicit Hamiltonian system as defined in section El Section El 
discusses the dynamics of a singular reduced implicit Hamiltonian system. It is shown that the 
"projectable" solutions of the original system project to solutions of the singular reduced system. 
Section El describes the decomposition of the singular reduced space into smooth manifolds called 
pieces. The singular reduced implicit Hamiltonian system restricts to regular reduced implicit 
Hamiltonian systems on the pieces. Section contains some concluding remarks. 
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2 Implicit Hamiltonian systems 



Let M be a smooth n-dimensional manifold and let TM © T*M denote the vector bundle whose 
fiber at x G M is T X M x T*M. (Throughout this paper all geometric objects are assumed to be 
smooth, so when manifolds, vector bundles, sections are mentioned, they are all smooth.) A Dirac 
structure on M is defined as follows. 

Definition 1. A smooth vector subbundle T> C TM ®T* M is called a Dirac structure if for every 
fiber V(x) C T X M x T*M, x G M, one has V(x) = V x (x), where 

V ± {x) = {(w,w*) G T X M x T*M | (v*,w) + (w*,v) = 0, V(v,v*) G T>(x)}. (2.1) 

(Here (•, •) denotes the duality pairing between TM and T*M .) 

Notice that T> being a vector subbundle of TM © T*M implies that its fibers all have the same 
dimension, i.e., dimD(x) = dimT>(x'), Vx,x' G M. In particular, if D is a Dirac structure then 
dim£>(x) = n, Vx G M. Furthermore, T>(x) = T> (x), x G M, implies that 

(v*,v) = 0, y(v,v*) G T>(x). (2.2) 

Remark 2. In |12j a constant Dirac structure on a vector space V is defined as a vector subspace 
V C V x V* such that V = V^. An equivalent way of writing Definition ^ is therefore: a Dirac 
structure on a manifold M is a smooth vector subbundle T> C TM © T*M such that each fiber 
T>(x), x G M, is a constant Dirac structure on T X M. ♦ 

In [B] there is yet another slightly different definition of a Dirac structure. Denote by %i oc (M) 
(respectively X(M)) the space of local (respectively global) smooth sections of TM, that is, these are 
the spaces of smooth local (respectively global) vector fields on M. Similarly, Q^ oc (M) and Q k (M) 
denote the spaces of smooth local and global fc-forms on M. The spaces of smooth local and 
global sections of the vector subbundle V C TM T*M are denoted by T>i oc and D respectively. 
Throughout, let X,Y G X ioc (M) and a, (3 G tt\ oc {M). Define a pairing on smooth sections of 
TM © T*M by 

(((X,a), {Y,/3))) = (a,Y) + {p,X), for(X,a), (Y, /3) G TM © T*M. (2.3) 
According to [H], a Dirac structure on M is a smooth vector subbundle T> C TM © T*M such that 

1. D is isotropic, i.e., for every two (local) sections (X, a), (Y, (3) G T>i oc : (((A, a), (Y, (3))} = 0; 

2. £> is maximal, i.e., if (Y,(3) is a (local) section of TM © T*M such that (((A, a), (Y,/3))) = 
0, V(A,a) G 2)* oc , then (Y, f3) G 2) /oc . 

It is easily shown that this definition given in |Sj and Definition Q are equivalent. Indeed, since T> 
is a smooth vector subbundle, every (v, v*) G T>(x) can be extended to a local section (A, a) G S); oc - 
Furthermore, T> being a smooth vector subbundle implies that also TJ^ C TM © T*M (with fibers 
P J -(x)) is a smooth vector subbundle and therefore also every (w,w*) G T> ± (x) can be extended 
to a local section (Y, (3) of 2? -1 -. Elementary linear algebra shows that 2? is isotropic if and only if 
T> C T) 1 - and that maximal isotropy is equivalent to equality in this inclusion or to the fact that 
dimD(x) = n for all x G M (see, e.g., pQ, Section 5.3 for this type of argument). 
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Corresponding to a Dirac structure a number of (co-)distributions may be identified. Recall 
that a distribution A on a manifold M is a map which assigns to each x G M a linear subspace 
A(x) of the tangent space T X M. A is called a smooth distribution if around any point these 
subspaces are spanned by a set of smooth vector fields, i.e., there exist X±, . . . ,Xk G Xi oc (M) 
such that A(x) = span {Xi(a;),... , -Xfc(x)}. The distribution A is called constant dimensional if 
the dimension of the linear subspace A (a;) C T X M does not depend on the point x G M. Notice 
that if A is a smooth constant dimensional distribution on M, then it defines a smooth vector 
subbundle (also denoted by A) of the tangent bundle TM, with fibers A (a;), x G M. Analogously, 
a codistribution T is defined as a map which assigns to each x G M a linear subspace of the 
cotangent space T*M. Smoothness and constant dimensionality are defined in the same way as for 
distributions. A smooth constant dimensional codistribution defines a smooth vector subbundle of 
the cotangent bundle T*M. 

Now, any Dirac structure T> naturally defines a distribution A, with fibers given by 1 

A{x) := {X(x) | X G Xi oc (M), (X, 0) G ® loc }, (2.4) 

and a codistribution T, whose fibers are defined by 2 

T(x) := {a(x) | a G nj oc (M), 3X G X loc (M), (X, a) G £ ioc }. (2.5) 

Since T> is isotropic it follows that A(x) C T°(x), where r°(x) denotes the annihilating vector 
subspace of T(x) in T X M, that is, T°(x) = {v G T^M | (v*,v) = 0, Vw* G r(a;)}. Equivalently, 
r(x) C A°(x), where A°(x) denotes the annihilating vector subspace of A(x) in T*M, that is, 
A°(x) := {v* G T*M \ (v*,v) = 0, Vw € A(x)}. Furthermore, if T is constant dimensional, and 
hence defines a vector subbundle of T*M, it follows by maximal isotropy of T> that A{x) = T°(x) 
(equivalently, T(x) = A°(x)). Notice that in this case also A is constant dimensional and hence 
defines a vector subbundle of TM. 

A special property of a Dirac structure is defined in the following 

Definition 3. A Dirac structure V is called closed, or integrable, if for all (local) sections (Xi,ai), 
(A 2 ,a 2 ), (A 3 ,a 3 ) G T>i oc 

(L Xl a 2 ,X 3 ) + (L X2 a 3 , X x ) + {L Xa a 1 ,X 2 ) = 0. (2.6) 
Equivalently T W^ . V is closed if and only if for all (Xx,a±), (X2,a 2 ) G Di oc 

([X 1 ,X 2 ],L Xl a 2 - Lx 2 ai + d(ax,X 2 )) G 2> ioc . (2.7) 

The notation L x is reserved for the Lie derivative operator (acting on any type of tensor field) 
defined by the (local) vector field X on M. 

It is easy to see that the graph of a symplectic form u : TM — > T*M or the graph of the skew- 
symmetric vector bundle map J : T*M — > TM induced by a Poisson bracket on M defines a Dirac 

1 In the literature on implicit Hamiltonian systems this distribution is usually denoted by Go and sometimes called 
the characteristic distribution. However, in order to avoid confusion with notation defined later on in this paper we 
decided to adopt a different notation here. 

2 This codistribution is usually denoted by Pi in the literature. 
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structure on M. As customary, we will call both the bundle map J and the two-tensor defined by 
{■, •} the Poisson structure on M. Closedness of T> corresponds in these two cases to the condition 
that a; is a closed two-form, respectively, the Poisson bracket satisfies the Jacobi identity. 

In this paper we will concentrate on a rather frequently occurring type of Dirac structure defined 
as follows. Let {•,•} : C°°(M) x C°°(M) -> C°°(M) be a generalized Poisson bracket on M. That 
is, {•, •} is skew-symmetric, bilinear, and satisfies the Leibniz property. Denote the corresponding 
vector bundle map by J : T*M -» TM (i.e., J(dH,dF) = {H, F}, H,F G C°°(M)); recall that J 
is skew-symmetric. Note that we do not require {•, ■} to satisfy the Jacobi identity and neither that 
J have constant rank. Moreover, given a subbundle A of TM (i.e., a smooth constant dimensional 
distribution A on M), it is easy to see that the vector subbundle T> C TM © T*M with fiber 

V{x) = {(v,v*) G T X M x T*M | v - J{x)v* G A(x), v* G A°(x)} (2.8) 

defines a Dirac structure on M. In terms of its local sections this is expressed as 

£>Zoc = {(X, a) G Xi oc (M) © O toc (M) | X — J a is a local section of A, a is local section of A°}, 

(2.9) 

where A° denotes the vector subbundle of T*M whose fiber at x G M equals A°(x). 

Remark 4. In ^5] it is shown that, under a mild constant dimensionality assumption, every Dirac 
structure can be written in the form (|2,8|) or, equivalently, Q2.9|) . Indeed, if T> is an arbitrary Dirac 
structure on M, define the codistribution T as in 1)2. 5 J) , and assume that T is constant dimensional 
and hence defines a vector subbundle of T*M. Then there exists a well defined (see |12| ) skew- 
symmetric vector bundle map J(x) : T(x) C T*M — > (r(x))* C T X M, x G M, defined by 

j( x ) v * = ve (r(x))*, v* g r(x), (2.10) 

where v G (T(x))* is the restriction of some v G T X M to T{x) C T*M which satisfies the condition 
that (v,v*) G V{x). Notice that the kernel of J{x) is given by the codistribution To with fibers 
defined by 3 

T (x) := {a(x) I a G nj oc {M), (0, a) G D /oc }. (2.11) 

Then the Dirac structure T> takes the form l|2.8|l . or, equivalently, (|2.9|) . with A = T°. 
The map J may be locally extended to a generalized Poisson structure J : T*M — ► TM. ♦ 

Notice that, in general, the Dirac structure defined in Q2.8JI is not closed. Although closedness 
does not play an important role in the rest of the paper we remark for completeness that (|2,9|) 
defines a closed Dirac structure if and only if (see |12j ) 

1. A is involutive, 

2. the bracket {•, •} restricted to the set of admissible functions Ad '■= {H G C°°(M) [ dH G A°} 
defines a Poisson bracket on Ad (that is, the Jacobi identity holds). 

Remark 5. Before leaving our introduction to Dirac structures and moving on to the description 
of implicit Hamiltonian systems, we would like to mention the following generalization. In JH] a 
3 This codistribution is usually denoted by Po in the literature. 
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Dirac structure is defined as a maximal isotropic subbundle V C A © A*, where the pair (A, A*) 
defines a Lie bialgebroid over a smooth manifold M . The isotropy is defined with respect to the 
natural pairing ((,)) defined analogously as in Q2.3|) by A and its dual A*. If we take the special 
case A = TM and dually .A* = T*M, then we are back at the definition given earlier in this paper. 
For more information on this generalization we refer to |15j and the references therein. We remark 
that 15 require the Dirac structure to be closed. In their terminology, Dirac structures are always 
closed, while Dirac structures not satisfying condition (|2.7|) are called almost Dirac structures. In 
this paper however, we prefer to use the terminology as introduced above, calling a maximally 
isotropic subbundle of TM © T*M a Dirac structure and adding the prefix closed if (and only if) 
the conditions of Definition |31 are satisfied. ♦ 

Now we turn to the definition of an implicit Hamiltonian system. Consider a Dirac structure T> 
on M and a smooth function H G C°°(M), called the Hamiltonian or energy function. Then the 
three-tuple (M, T>, H) defines an implicit Hamiltonian system as follows: 

Definition 6. The implicit Hamiltonian system (M, T>, H) is defined as a set of smooth time 
functions {x{t) \ x : M — > M of class C°°} (called solutions ) satisfying the condition 

(x(t),dH(x(t))) G V(x(t)), Vt. (2.12) 

Equations (|2.2JI and ([2.12ft imply that implicit Hamiltonian systems are energy conserving, i.e., 

dH 

— (x(t)) = (dH(x(t)),x(t)} = 0, Vt (2.13) 

If T> is the graph of a symplectic form a; or of a Poisson structure J : T*M — > TM then the above 
definition becomes that of a classical (explicit) symplectic or Poisson Hamiltonian system. On the 
other hand, if V is defined by (|2.8[) then the system includes the algebraic constraints 

dH{x{t)) G A°(x(t)), Vt. (2.14) 

Thus all solutions of the implicit Hamiltonian system necessarily lie in the constraint manifold 

M c = {xeM\ dH{x) G A°(x)}. (2.15) 

Since the implicit Hamiltonian system defines a set of differential and algebraic equations, there 
is not an existence and uniqueness result as one has for solutions of classical Hamiltonian systems 
described by ordinary differential equations. In particular, not every point xq G M c necessarily lies 
on the trajectory of some solution x(t) of the system, and neither are the solutions through a point 
xq G M c (if existing) necessarily unique (this happens, for instance, if the Lagrange multipliers 
corresponding to the algebraic constraints cannot be solved uniquely). In the sequel we will not 
investigate these problems. Instead we will study the reduction of the underlying Dirac structure 
in the presence of symmetries (defined later on) and show that certain "projectable" solutions (if 
existing) will project to solutions of an implicit Hamiltonian system on the reduced space. 

The problem of existence and uniqueness of solutions to implicit systems is an important and 
active area of research and will not be touched upon here. We would only like to mention the 
special case of so-called index 1 systems. Consider the implicit Hamiltonian system defined by 
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the Dirac structure ()2.8I) and the Hamiltonian function H € C°°(M). Let the vector subbundle A 
be (locally) written as the span of the independent vector fields g± , . . . , g m . Then the constraint 
manifold can be written as 



M c = {x e M | L g H{x) = 0, j = 1, . . . , m}. 



(2.16) 



Now assume that the constraints are of index 1, that is, the matrix 



(2.17) 



is nonsingular for all x € M c . In that case one can restrict the implicit Hamiltonian system 
(M, T>, H) to an explicit Hamiltonian system on M c , defined by a (possibly nonintegrable) Poisson 
bracket on M c , see |32| I5]. Its corresponding dynamics is thus given by a set of ordinary differential 
equations on M c . Standard existence and uniqueness results now yield that through every point 
xq G M c there is a unique solution of the implicit Hamiltonian system (restricted to M c ). 

Example 7. As an important example of implicit Hamiltonian systems we like to mention the 
class of mechanical systems with kinematic constraints. These systems are described by implicit 
Hamiltonian systems (M, T>, H) with T> of the form (|2,8|) . Here, H is the total energy of the system, 
(the sum of kinetic and potential energies for classical mechanical systems) and the phase space 
M = T*Q is the cotangent bundle of the configuration space Q; local coordinates are denoted, as 
usual, by (q,p) € T*Q. The Poisson bracket {•,•} is the standard Poisson bracket corresponding 
to the canonical symplectic form u = dp A dq on T*Q (i.e., the associated Poisson structure J 
is given by lo^ 1 ). Finally, if we assume the kinematic constraints to be linear in the velocities 
(e.g., non-slipping constraints of a rolling wheel) then they can be (locally) written in the form 



By d'Alembert's principle, the constraints (|2. 18|) generate constraint forces of the form F c = A(q)X, 
where A are the Lagrange multipliers. It follows that the distribution A is (locally) described by 
the image of the matrix 



If the kinetic energy is defined by a positive definite metric on Q, then the constraints are of index 
1, i.e., the Lagrange multipliers A can be solved uniquely. In this case the constrained mechanical 
system on T*Q can be written as an unconstrained generalized Hamiltonian system on M c . In 30 
it is shown that the corresponding Poisson bracket on M c satisfies the Jacobi identity if and only 
if the kinematic constraints are holonomic. 

3 Symmetries and regular reduction 

In this section we recall some of the results in [SJ El concerning symmetries and reduction of implicit 
Hamiltonian systems. These results act as a reference for the results obtained in this paper and 
will be specifically needed in section H3 to show that the singular reduced Dirac structure restricts 
to regular reduced Dirac structures on the pieces corresponding to the orbit type decomposition 



A T {q)q = A T 




(2.18) 
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of the singular reduced space Mq . We refer to EJ for a detailed treatment of symmetries and 
reduction of implicit Hamiltonian systems. We stress that, unless specifically stated otherwise, the 
results in this section are valid for arbitrary Dirac structures, not necessarily of the form (|2.8[) . 

Definition 8. A smooth vector field Y on M is called a symmetry of the Dirac structure T> if for 
every (local) section (X, a) ofT>, one also has that (LyX, Lya) is a (local) section ofD. Y is 
called a symmetry of the implicit Hamiltonian system (M,T>,H) if Y is a symmetry of T> and a 
symmetry of H, i.e., LyH = 0. 4 

The above notion of symmetry generalizes the classical notion of symmetry in symplectic or 
Poisson Hamiltonian systems. Indeed, if T> is the graph of a symplectic form u then Y is a symmetry 
of T> if and only if Lyuj = 0. Likewise, if T> is the graph of the skew-symmetric vector bundle map J 
defining a Poisson structure on M, then Y is a symmetry of T> if and only if the integral flow of Y is a 
Poisson map, i.e., Y is a derivation of the Poisson bracket: Ly{F\, F2} = {LyFi, F2} + {F±, Lyi 7 ^} 
for any F 1 ,F 2 £ C°°(M). 

In the next sections we will consider a special subclass of symmetries defined by the following 

Proposition 9. Consider a Dirac structure T> of the type defined in \2. 8\) . If the vector field Y on 
M is a derivation of the generalized Poisson bracket (equivalently, LyJ = 0) and LyZ is a local 
section of A whenever Z is a local section of A, then Y is a symmetry ofT>. 

In particular this means that we restrict our attention to the case where Y is a symmetry of the 
generalized Poisson bracket as well as a symmetry of the vector subbundle A C TM. These kinds 
of symmetries often arise in constrained mechanical systems, see also Remark |^] later on. 

Specifically we will consider Lie algebra symmetries, defined as follows. Recall that a (left) Lie 
algebra action on a manifold M is a Lie algebra anti-homomorphism £ £ q i— * £m £ 3£(M) such 
that the map (x,£) £ M xjh £m( x ) £ T X M is smooth. Then q is a symmetry Lie algebra of 
T> if £jvf is a symmetry of T> for every £ £ Q. In particular, if the Dirac structure is of the type 
defined in (|2.8[) . the criterion in Proposition applies. Similarly, q is a symmetry Lie algebra of an 
implicit Hamiltonian system (M, T>, H) if each £m is a symmetry of T> and also a symmetry of H, 
i.e., L^ M H = 0. 

Lie algebra symmetries are often induced by Lie group actions. If G is a Lie group with Lie 
algebra g and (j) : G X M —> M is a smooth left action of G on the manifold M, the infinitesimal 
generator of the action associated to £ £ g defined by 

Zm(x) = j t <P(eM&),x)\t=o £ T X M, x £ M, (3.1) 

induces a left Lie algebra action of q on M. Then G is said to be a symmetry Lie group of T> if q 
is a symmetry Lie algebra of T>. Similarly, G is a symmetry Lie group of the implicit Hamiltonian 
system (M, T>, H) if q is a symmetry Lie algebra of this implicit Hamiltonian system. 

We turn now to the analysis of the regular reduction process of Dirac structures and implicit 
Hamiltonian systems. We start by explaining how an implicit Hamiltonian system on M can be 
restricted to an implicit Hamiltonian system on a submanifold N of M. Let T> be & Dirac structure 
on M and let N C M be a submanifold of M. Following 8 define for each x £ N the map 

4 This is called a strong symmetry in 
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a(x) : T X N x T£M -» T X N x T x *iV, x G JV, by a(a?) («,«*) = (v, v* \ TxN ), where v*|t,tv means the 
restriction of the covector v* G T*M to the subspace T X N C T X M. Define a vector subspace of 
T X N x T X *JV by 

2^(3) = a(x) (V(x) n (T^JV x T*M)) , x £ N. (3.2) 

It is clear that T>n(x) C T>j^(x), x G iV. To prove the reverse inclusion, suppose that (w,w*) G 
£>^(x) C T^JV x T*JV, i.e., («*,«;) + («>*,«) = 0, V(«,«*) G X»jv(x). Then = a(x)(v,v*), 

where v* G T*M, (v,v*) G 2?(x), ^It^n = v* and since v,w £ T X N, one gets 

= (v*,w) + (w*,v) = (v*,w) + (w*,v), 

where w* is an arbitrary extension of w* to T X M. Since this relation holds for all (v,v*) G T>{x) 
with v G T^JV, this implies that 

(it;, tD*) G [V{x) n (T X N x r*M)] x = V L {x) + (T^JV x t'm)- 1 

= V(x) + ({0} x T X N°), (3.3) 

so there exists a u* G T X N° C T*M such that (w,w* + u*) G V{x). However, since w G T X N and 
<t(x)(«;, + u*) = (w, (w* + u*)\t x n) = {w, w*), it follows that (w,w*) G T>n{x), which shows that 
£> jv(sc) C T>n(x). Assuming that the dimension of V(x) D (TajJV x T*M) is independent of x G JV, 
that is, that £> n (TN x T*M) is a vector subbundle of TJV x T*M, it follows that a is a vector 
bundle map and hence that T>n is a vector subbundle of TN x T*N. So we have proved (a slightly 
rewritten version of jSJ): 

Proposition 10. Consider a Dirac structure T> on M and let N be a submanifold of M . Assume 
that T>{x) n (T X N x T*M), x G JV, has constant dimension on JV. Then the bundle T>n with fibers 
defined by is a Dirac structure on N . This is called the restriction ofT> to JV. 

In order to do computations it is convenient to describe the restricted Dirac structure T>n in 
terms of its local sections. This gives the following proposition (an improved version of [31|S]). Let 
1 : JV <^-> M denote the inclusion map. 

Proposition 11. Consider a Dirac structure D on M and let N be a submanifold of M. Assume 
that T>{x) H (T X N x T*M), x G JV, has constant dimension on JV and let T>n denote the restriction 
ofT> to JV. Then (X,a) is a local section ofD^ if and only if there exists a local section (X, a) of 
T> such that X ~ t X and a = L*a. Otherwise stated, in terms of its local sections 

(£> n)ioc = {(X, a) G Xioc(N) ^l oc (N) I 3(X, a) G Tl loc such that X ~ t X and a = L*a}. (3.4) 

Here X ~ t X denotes the fact that X and X are t-related and 1* denotes the pullback by 1. It 
can be shown that if T> is closed, then also D/v is closed. 

Now, let (M, T>, H) be an implicit Hamiltonian system on M and let JV be a submanifold of 
M such that the constant dimension condition of Proposition EH is satisfied and assume that (the 
flow corresponding to) the solutions of (M,T>,H) leave the submanifold N invariant. Restrict the 
Hamiltonian H to a smooth function on N, i.e., Hn = Hol, and define the implicit Hamiltonian 
system (N,T>n, Hn) on JV. Then: 
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Proposition 12. Every solution x(t) of (M,D,H) which is contained in N is a solution of 
(N,V N ,H N ). 

We remark that, in general, there is not a one-to-one correspondence between the solutions gen- 
erated by the original system (M, V, H) and those generated by the restricted system (N, T>n, Hjy). 
(Indeed, compare this with the case of restriction of a symplectic form uj on M to an arbitrary 
submanifold N, leading to a nontrivial kernel for ujjy.) In case N happens to be the level set of a 
Casimir function of the Dirac structure T> then there is a one-to-one correspondence between the 
solutions of the original system and of the restricted system; see [3] for more information. 

Next, we explain how an implicit Hamiltonian system on M admitting a symmetry Lie group 
G can be projected to an implicit Hamiltonian system on the orbit space M/G. Consider a Dirac 
structure T> on M and let G be a symmetry Lie group of T>, acting regularly on M, that is, 
the orbit space M/G is a smooth manifold and the canonical projection map tt : M — > M/G 
is a surjective submersion. Let V = kerT-7r denote the vertical subbundle of TM, with fiber 
V{x) = span {£m(x) I £ G fl} for every x G M. We assume that V + A is a smooth vector 
subbundle of TM, i.e., its fibers all have the same dimension. Furthermore, define the smooth 
vector subbundle E C TM © T*M in terms of its local sections by 

T loc (E) = {(X, a) G Xi oc (M) © Qj oc {M) | a = it* a for some a G ttj oc (M/G)} (3.5) 

(where Ti oc (E) is the space of local sections of the subbundle E) and assume that T> HE is a smooth 
vector subbundle of TM © T*M, i.e., its fibers all have the same dimension. Then we have: 

Proposition 13. 1281 0/ Consider a Dirac structure V on M admitting a symmetry Lie group G 
acting regularly on M . Assume that V + A is a smooth vector subbundle ofTM and that Vd E is 
a smooth vector subbundle of TM ®T* M . Then T> "projects to a Dirac structure T> on M := M/G, 
described in terms of its local sections by 

®ioc = a) £ Xi oc (M) x 9} loc (M) | 3(X, a) £ D loc such that X X and a = ir*a}. (3.6) 

T> is called the projection of T> to M/G. Again, it can be shown that closedness of T> implies 
closedness of V. Let (M, T>, H) be an implicit Hamiltonian system admitting a symmetry Lie group 
acting regularly on M such that the conditions in Proposition are satisfied. The G-invariant 
function H defines a function H G C co (M/G) by H = H o tt. Consider the implicit Hamiltonian 
system (M/G,T>,H). A G-projectable solution x(t) of (M,T>,H) is defined as a solution x(t) of 
(M,V,H) for which there exists a projectable vector field X G Xi oc (M) (i.e., X ~„- X for some 
X G 3Li oc (M/G)) such that x(t) = X(x(t)). The following proposition was obtained in 

Proposition 14. If x{t) is a G-projectable solution of (M,T>, H) then x(t) := n(x(t)) is a solution 
of (M/G,T>,H). Conversely, every solution x(t) of (M/G,T>,H) is locally the projection under n 
of a G-projectable solution x(t) of (M, T>, H). 

In [5] a simple example is given showing that not every solution of an implicit Hamiltonian 
system (M, T>, H) admitting a symmetry Lie group G is necessarily G-projectable. However, if 
the constraints are of index 1, and therefore the implicit Hamiltonian system can be restricted to 
an explicit Hamiltonian system on the constraint manifold M c , then it can be shown that every 
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solution is G-projectable and therefore projects to a solution of the reduced implicit Hamiltonian 
system on M/G. 

Finally, let us briefly recall some results in [SJ El concerning the reduction of implicit Hamiltonian 
systems admitting a symmetry Lie group having an associated momentum map. Consider an 
implicit Hamiltonian system (M, T>, H) admitting a symmetry Lie group G with Lie algebra 0. 
Assume that there exists an Ad* -equi variant map P : M — > g* (Ad* denoting the coadjoint action), 
called momentum map, such that 

(ZM,dPt) G2, V£G , (3.7) 

where P$ G C°°(M) is defined by P^x) = (P(x),£), x G M. Notice that if V is the graph of a 
symplectic form w or a Poisson structure J, then this corresponds to the classical definition of a 
momentum map. Assuming that fi G g* is a regular value of P, it follows that the level set P _1 (/x) 
is a closed submanifold of M. Since the Hamiltonian is G-invariant, the solutions of (M, T>, H) leave 
the level set P _1 (/x) invariant. Thus by ()2.12j) . (|3.7j) . and the identity V(x) = T> ± (x) it follows that 

^(x(t)) = (dPt(x(t)),x(t)) = -(dH,^ M )(x(t)) = 0, Vt, V£ G 0. (3.8) 

In other words, P is a first integral of the implicit Hamiltonian system (M,T*,H). Assuming that 
the conditions in Proposition ^] hold, we can restrict the implicit Hamiltonian system (M, T>, H) 
to an implicit Hamiltonian system (N,T>n, Hn) on N = P^ 1 (fi). The system (N,T>n, Hn) admits 
the symmetry Lie group G M := {g G G \ Ad*^ = fx}. Assume that G^ acts regularly on N, that is, 
N/G^ is a smooth manifold with the canonical projection a surjective submersion, and assume that 
the conditions in PropositionEDare satisfied. Then we can project the implicit Hamiltonian system 
(N,T>n, Hn) to an implicit Hamiltonian system (M^, P M , H^), where M M = N/G^ = P~ 1 (fj,)/G fl is 
the regular reduced space and G C°°(M M ) defined by on = Hn is the reduced Hamiltonian. 
The implicit Hamiltonian system (Mu,'Du>Hu,) is called the reduced implicit Hamiltonian system 
corresponding to (M, V, H) and T>^ is called the reduced Dirac structure. If V is the graph of a 
symplectic form u, then is precisely the graph of the Marsden-Weinstein reduced symplectic 
form Uu. Likewise, if T> is the graph of a Poisson structure J on M, then is the graph of the 
reduced Poisson structure |17j . Notice, however, that contrary to the above mentioned classical 
reduction results, closedness of the Dirac structure is not required for the reduction scheme to 
work. This observation is important since it allows the reduction scheme to be applied to the class 
of mechanical systems with (possibly nonholonomic) kinematic constraints; see [5J |S] for further 
information and a discussion of its relationship to other recent results in this area. Of course, if T> 
happens to be closed then also the reduced Dirac structure T>^ will be closed. 

Finally, we would like to mention that in [SJ El it is shown that the reduction scheme can also 
be applied the other way around, starting with factorizing the symmetry group G and afterwards 
restricting the result to a level set of the remaining first integrals (which actually turn out to be 
Casimir functions). The resulting implicit Hamiltonian system on is isomorphic to the system 
(Mu,Pu, Hu)- Notice that is actually the orbit reduced space P- 1 (0^)/G, where C M denotes 
the coadjoint orbit in 0* through fi. See also El 1221 for some classical references. 

Intrinsic reductions In this paragraph we will set the reduction results described above against 
what we will call intrinsic reductions. The latter are independent of any symmetry properties of 
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the Dirac structure and, in fact, can be perfomed on any closed Dirac structure. These kind of 
reductions have been described in the literature by various authors [HI El El El El • 

Consider a closed Dirac structure PonM. Then by condition (|2.7|) it follows that the charac- 
teristic distribution A is involutive, and hence by Frobenius' Theorem defines a regular foliation 
of M into integral submanifolds of A. On the other hand, the distribution defined by 5 

9(x) := {X{x) | X G X loc (M), 3a G nj oc (M), (X, a) G V loc } (3.9) 

clearly is also involutive, defining a regular foliation <l?e of M into integral submanifolds of B. (For 
the moment we will assume that both distributions are constant dimensional.) 

There a two logical ways to "reduce" the Dirac structure on M to a lower dimensional manifold. 
The first is to project the Dirac structure to the quotient manifold M/<3?a ; he., by factoring out 
the characteristic distribution. This was done in (Sj where it was shown that the Dirac structure T> 
on M induces a well defined Poisson bracket on the quotient manifold M/<I>a (0; Corollary 2.6.3). 
This remarkable result was generalized in to Dirac structures on Lie bialgebroids as described 
in Remark El where is was refered to as Poisson reduction. In 5 j it was observed that in fact this 
reduction can be seen as a special case of symmetry reduction if one notices that the distribution A 
is a symmetry distribution of T>, i.e., every vector field Y G A is a symmetry of T> as in Definition 
El The Dirac structure T> can be projected to a Dirac structure T> on M/<I>a using Proposition 
1131 It turns out that T> is exactly the graph of the Poisson structure corresponding to the Poisson 
bracket defined by Courant |S]. We refer to [5], Example 4.2.4, p. 73, for more details. 

The second reduction possibility is to restrict the Dirac structure to each of the integral subman- 
ifolds of B. This can be done using Proposition El and results in a Dirac structure on each of the 
integral submanifolds of 6. In 0E] it is shown that each of the reduced Dirac structures represents 
a presymplectic structure on the corresponding leaf of the foliation, see p], Example 4.1.8, p. 69, 
and [6_, Example 9, p. 79. This corresponds to Theorem 2.3.6 in Courant [H] and Theorem 2.2 in 
Dorfman stating that a closed Dirac structure has a foliation by presymplectic leaves. 

Once more we want to stress that the reductions described above are "intrinsic" and have nothing 
to do with the existence of any symmetry groups of the implicit Hamiltonian system (although, 
as explained above, the first reduction can be interpreted in terms of symmetries of the Dirac 
structure). They can be perfomed on any closed Dirac structure. We will not concentrate on 
these intrinsic reductions anymore, and instead will investigate symmetry Lie groups of implicit 
Hamiltonian systems, together with their (singular) reductions. Doing so, we do not assume that 
the Dirac structure is closed and in fact all our results will be presented for the general case. 

4 Singular reduction 

Contrary to the regular reduction reviewed in the previous section we now describe in a purely 
topological way how to obtain a reduced Dirac structure on the reduced space M M if M M is not a 
manifold. This is the case when /x is a singular value of the momentum map P. In that case, vector 
fields and differential one-forms on are not defined and therefore the results described in the 
previous section cannot be used. Describing the dynamics corresponding to such a topologically 

5 This distribution is usually denoted by Gi in the literature. 
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reduced Dirac structure on M^ will be done in section |SJ For easiness of exposition we will take 
\x = throughout the rest of this paper. 

From this point on we specifically consider only Dirac structures of the form (|2.8j) . admitting 
symmetries as described in Proposition |5J that is, given a vector subbundle A C TM, and a 
generalized Poisson structure J : T*M — > TM (the Jacobi identity does not necessarily hold), the 
Dirac structure is defined by 

V(x) = {(v,v*) G T X M x T*M | v - J(x)v* G A(x), v* G A°(x)} (4.1) 

and Y G 3C(M) is a symmetry of T> if 

LyJ = and LyZ is a local section of A whenever Z is a local section of A. (4.2) 

Consider such a Dirac structure P on a manifold M admitting a symmetry Lie group G with 
corresponding A<i*-equivariant momentum map P satisfying ([3.7)1 : it is not assumed that G acts 
regularly on M. Let \i = G 0* be a singular value of P and consider the level set N = P _1 (0) 
which is not a smooth submanifold of M. However, TV is a closed subset of M and is a topological 
space relative to the induced subspace topology. The level set N is G-invariant so one can endow 
the orbit space Mo := N/G = P _1 (0)/G with the quotient topology. Denote by tt : N — > Ma the 
canonical projection map, that is, tt maps x G N onto its orbit G ■ x G Mq. 

Define the set of smooth functions on Mo as follows (see [0], or for the original source see [25] ) . 

Definition 15. A continuous function /o on Mq is called smooth, denoted by /o G C°°(Mo), if 
there exists a smooth G-invariant function f G C OD (M) G such that fo°^ = /|p-i(o)- 

Given the singular reduced space Mo together with its topology and a set of smooth functions 
C°°(Mo) on Mo, we want to define a reduced Dirac structure on Mo. Let P be a Dirac structure 
on M of the type defined in ([4.1)1 and assume that the infinitesimal generators £ G Q, satisfy 
the conditions in ([4.2)1 . Since G is a symmetry Lie group of the generalized Poisson bracket {•, •} : 
C°°(M) x C°°(M) — > C°°(M), corresponding to the bundle map J, we can use the theory in jHUH] 
to define a generalized Poisson bracket {-,-}o : C°°(Mq) x C°°(Mo) — > C°°(Mo) on the singular 
reduced space M . This goes as follows. Let f ,h G C°°(Mo) and let f,h £ C°°(M) G be such 
that /o o tt = /|p-i(o) and ho o tt = h\p-ir \. Define the singular reduced generalized bracket by 

{/o,Mo 07r = {fih}\ P -i(py (4.3) 

This gives a well defined generalized Poisson bracket on Mo. In particular, ([4.3)1 does not depend 
on the choice of the G-invariant extensions / and h (whose existence is assumed, by definition). 

Remark 16. The reduction theory in [HJIH] is only developed for the singular reduction of symplectic 
manifolds under a symmetry Lie group action. That is, the Poisson bracket {•, •} is assumed to 
be nondegenerate and to satisfy the Jacobi identity. In principle, however, these results generalize 
immediately to the case of singular reduction of generalized Poisson brackets, as described above by 
([4.3)1 . In particular, [3110] show that (under the assumption that G acts properly) nondegeneracy of 
{•, •} implies that of {•, -}o- Also, from ([4.3)1 it follows immediately that {•, -}o satisfies the Jacobi 
identity if {•, •} does. 

Once again for clarity: In this paper we do neither assume that the generalized Poisson bracket 
{•, •} is nondegenerate, nor that it satisfies the Jacobi identity. Furthermore, properness of the 
group action is not assumed until Section IFJJ ♦ 
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Next, consider the vector subbundle A C TM, denning a constant dimensional distribution on 
M. We show that A defines a vector space A consisting of derivations on the space C°°(Mq) of 
smooth function on Mq. Denote by r; oc (A) the local sections of the subbundle A C TM. We show 
that every vector field X G Ti oc (A) is "tangent" to N. In the regular case, when N = is 
a smooth submanifold of M, this means that X restricts to a well defined vector field X on N. 
However, if = is a singular value of the momentum map, then N is not a smooth manifold 
and hence we have to define what "tangent" means. Recall that a vector field X G X{M) is in 
one-to-one correspondence with a derivation, also denoted by X : C°°(M) — * C°°(M), on the set 
of smooth functions on M. The correspondence is given by the formula 6 

X[f] = (df,X), V/gC°°(M). (4.4) 

A derivation X on M is said to be tangent to N if it restricts to a well defined derivation X on 
the set of Whitney smooth functions on iV. A continuous function / on N is said to be a Whitney 
smooth function if there exists a smooth function / on M such that / = /|jv; the set of Whitney 
smooth functions on N is denoted by W°°(N). Otherwise stated, X is tangent to N if there exists 
a derivation X on W^iN) such that X[f](x) = X[f\ N ](x), Vx G N, for all f G C°°{M). A 
necessary and sufficient condition for X to be tangent to N is that 

X[f](x) = X[h](x), VxGiV, (4.5) 

for all f,h€ C°°{M) such that /|jv = ^|jv- Notice that in case N is a smooth submanifold of M 
and N is closed in M, then the set W°°(N) of Whitney smooth functions on N is equal to the set 
C°°(N) of all smooth functions on N (as defined by the differential structure on the submanifold 
N), and the above given definition yields the usual meaning of a vector field being tangent to the 
submanifold N (and consequently restricting to a vector field on N). 

Consider a vector field (or equivalently, derivation) X on M, and define j(t) to be an integral 
curve of X through xq G M if 7 

j t f{l{t)) = X[/]( 7 (t)), Vt, V/ G C°°(M), 7 (0) = x . (4.6) 
Now, let X G r/ oc (A) and j(t) be an integral curve of X through xq G P _1 (0). In particular, 

j t Pdl(t)) = X[P 6 ](j(t)) = 0, Vt, V£ G 0, (4.7) 

since by (|2.8|) (or 1)4.1(0 and (|3.7|) . dP^{x) G A°(x), Vx G M. This implies that the integral curve 
of X G r/ oc (A) through every xq G P~ 1 (0) is contained in P _1 (0) (conservation of the momentum 
map). Then by the equivalence of derivations and velocity vectors (remember that M is a smooth 
manifold) it follows that 

X[f](x ) = J t fm)\t=o = jh{l(t))\t=o = X[h](x ), (4.8) 

for all f,hE C°°(M) such that /|jy = h\]y. So we have shown that every vector field X G r; oc (A) 
is tangent to iV = P _1 (0) and consequently restricts to a well defined derivation X on W°°(N). 

6 The derivation is usually called Lie derivative and is also denoted by X[f] = Lxj '■ 
7 Take the coordinate functions f = x* to obtain the usual definition 7 l (t) = X l (-y{t)). 
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In conclusion, the constant dimensional distribution A on M restricts to a vector space A of 
derivations on W°°(N). If A is locally spanned by the independent vector fields X\, . . . ,X m , then 
A is locally spanned by the independent derivations X±,... ,X m . 

Using the above results we are able to show that the distribution A on M projects to a well 
defined vector space A of derivations on the smooth functions C°°(Mq). A vector field X on M is 
said to project to Mq if there exists a derivation X on C°°{Mq) such that for every / G C°°(M) , 
= X[fo](iv(x)), Vx G N, with fo defined by fo o ir = /|jy- It is clear that X restricts to 
a well defined derivation X on C°°(Mq) if and only if X[f](x) does not depend on the extension 
of fo o ir off iV to M and furthermore X[f]{x) = X[f]{y) for all x,y G N such that ir(x) = n(y). 
Now let X be a local section of A. Since X is tangent to N it follows that X[f](x) = X[f\jy](x) = 
X[fo o 7r](x), Vx G iV, and therefore its value does not depend on the extension of fo o ir off N to 
M. It remains to show that 

X[f](x) = X[f](y), Vx, y & N such that vr(x) = Tt(y). (4.9) 

In general, this is not true for every local section X of A. However, we will show that there exists 
a basis of local sections X\, . . . , X m of A which satisfies (j4.9|) . 

Denote by 2J; oc the space of local sections of the vertical distribution V defined by V(x) := 
span {£,m(x) | ^ € g}. Since Lg A/ r/ oc (A) C Ti oc (A) for every £ G g by (|4.2|) . it follows that 
[r; oc (A), 2J; oc ] C 53/ oc + r; oc (A). Indeed, taking an arbitrary local section of V of the form Y = 
^2ihi^ M , hi G C°°(M), i = l,...,r = dim g, £m,--->£m being a local basis of V, and letting 
X G r ioc (A), it follows that 

r r 

[X, Y] = [X, Y, hiCu] = E + (^x^)d/ e r ioc (A) + 2J /oc 

i=l i=l 

which proves the inclusion [r toc (A), 2J/ oc ] C 5J; oc + r ;oc (A). 

Assuming that the distribution V + A has constant dimension on M, the above inclusion implies 
that there exists a basis X\, . . . , X m of local sections of A such that [Xi, 2J/ c] C %3ioc, i = 1, ■ ■ ■ ,m; 
see e.g. Theorem 7.5 on page 214 in the book by Nijmeijer and van der Schaft |21| (in the notation 
of that theorem: the involutive distribution D is V, the distribution G is A, and one takes / = 0). 
In particular, [JQ,£m] G 2Jz c; which implies that for all / G C°°(M) G 

= [Xi,£ M ][f\ = Xi[L iM f] - L^{Xi[f\) = -L iM (Xi[f}), V£ G g. (4.10) 

This means that the function Xi[f] is G-invariant and therefore satisfies (|4.9|) . In conclusion, 
there exists a basis X\, . . . , X m of local sections of A such that each Xi projects to a well defined 
derivation Xi on C°°(Mo). The derivations X\, . . . ,X m locally span (in other words, form a basis 
of) a vector space of derivations on C°°(Mo), denoted by A. 

Remark 17. In the regular case, i.e., when fi = is a regular value of the momentum map and 
G acts freely and properly on M, the reduced space Mo is a smooth manifold. Furthermore, the 
set of smooth functions C°°(Mo) equals the set of smooth functions as defined by the differential 
structure on Mq. Indeed, since N = P~ 1 (Q) is closed in M, every smooth G-invariant function 
on N can be smoothly extended to a G-invariant function on M (Hj. In that case the notion of a 
"projecting derivation" as defined above has the usual meaning of projection of a vector field on 
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M to a vector field on the reduced space Mq. In particular there exists a basis X\, . . . , X m of local 
sections of A, which are tangent to N, such that the restrictions X±, . . . ,X m project to Mq. That 
is, each X{ is 7r-related to a vector field Xi on Mq. The projected vector fields X\, . . . , X m form a 
basis of local sections of A. ♦ 

So far we have defined a set of smooth functions C°°(Mq) on the singular reduced space Mq, 
together with a generalized Poisson bracket {•, -}o and a vector space A of derivations on C°°(Mq). 
Recall that the original Dirac structure T> (of the type given by (|2.8I) or (|4.1|0 on the manifold M 
was completely determined by the generalized Poisson bracket {•,•}, corresponding to J, and the 
distribution A. Therefore it makes sense to define a reduced Dirac structure on Mq as follows: 

Definition 18. Consider the singular reduced space Mq together with the set of smooth functions 
C°°(Mq). The singular reduced Dirac structure Vq is defined as the pair ({•, -}q, A). 

We also call T>q a topological Dirac structure. It will be shown in the next section that the singular 
reduced Dirac structure T>q defines a Hamiltonian formalism on the singular reduced space Mq. 

In order to better comprehend the meaning of the singular reduced Dirac structure defined in 
Definition we show that in the case of regular reduction the topological Dirac structure T>q 
exactly defines the regular reduced Dirac structure on Mq. 

Regular reduction. Suppose that = is a regular value of the momentum map and G acts 
regularly on M (that is, M/G is a smooth manifold and the projection M —* M/G is a surjective 
submersion; for example, if G acts freely and properly, these conditions are satisfied). According 
to the results described in Section |3] the Dirac structure T> on M is reduced to a Dirac structure V 
on the manifold Mq in two steps: firstly, T> is restricted to a Dirac structure T> n on N = P _1 (0) 
defined by (|3.4|) and, secondly, Vjy is projected to a Dirac structure V on Mq defined by (|3.6|) (with 
T>i oc replaced by (T>n)ioc)- Otherwise stated, in terms of its local sections, 

i) loc = {(X, a) E Xi oc (Mo) nj oc (M ) I 3(X, a) £ £ ioc such that X is tangent to N 

and l*X ~ 7r X, l*ol = 7r*d}, (4-11) 

where i : N M is the inclusion and l*X denotes the push forward of X to a vector field on N 
(that is, l*X is simply the restriction of X to N which is possible since N is a closed submanifold 
of M and X is tangent to iV by hypothesis). 

Consider the topological Dirac structure T>q given by Definition 1181 Since Mq is a manifold, the 
generalized Poisson bracket {•, -}o on the set of smooth functions C°°{Mq) (see also Remark 1 17(1 de- 
fines a skew-symmetric vector bundle map Jo : T*Mq — > TMq by Jo(dfo, dho) = {/o, /io}o> fo, ho € 
C°°(Mq). The vector space A of derivations on C°°(Mo) defines a constant dimensional distribu- 
tion of vector fields on Mq (in other words, a vector subbundle of TMq), also denoted by A. Then 
the topological Dirac structure T>q defines a Dirac structure on the manifold Mq, also denoted by 
T>q , which in terms of its local sections is given by 

CDoW = {(X, a) G X loc (M ) nj oc (M ) \X-J ae r /oc (A), a G T loc (A )}. (4.12) 

Indeed, T>q is a Dirac structure on Mq as defined in Definition ^ (notice that it is of the same form 
as in H2.9|) ). We show that T> = T>q. Since both are Dirac structures and therefore their fibers are 
of the same dimension (i.e., dim Mq), it is enough to show that T> C T>q. 
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If (X, a) is a local section of T>, then there exists a local section (X, a) of T> such that l*X r^ n X 
and t*a = 7r*a. Since (X, a) is a local section of T> one has 

Z := X — Ja is a local section of A, a is a local section of A°. (4-13) 

Consider the vector field Ja G X; oc (M). Since (Ja,a) G T>i oc it follows from ()3.7j) and D = X? 1 - 
that 

(Ja)[P^(x) = (dPt,Ja)(x) = -{a,£ M )(x) = -(a, 0)(tt(x)) = 0, Vx G JV, V£ G 0. (4.14) 

This implies that the vector field Ja is tangent to N. Furthermore, by construction of the reduced 
generalized bracket ()4.3|) it follows that L*(Ja) Jq6. Since also l*X ^ w X, equation (|4.13f) 
implies that there exists a vector field Z G Xi oc (Mq) such that l*Z ^ n Z. It follows that Z G r; oc (A) 
by construction of A. This yields 

X-J a = Z£ T loc (A). (4.15) 

By construction, the distribution A is spanned by vector fields Zi, . . . , Z m for which there exists a 
basis of vector fields Z±, . . . , Z m G r; oc (A) such that l*Zj ~ 7r Zj, j = 1, . . . , m. Since i*a = -K*a 
and a G r; oc (A°), it follows immediately that 

(6c, Zj} o 7T = (a, Zj) o i = 0, j = 1,..., m, (4.16) 

and therefore a G r; oc (A°). In conclusion, (X, a) is a local section of T)q. So we have shown that 
T> C T>q and since both are Dirac structures on Mq this implies that T> = T>q. 

We conclude that in the case of regular reduction the topological Dirac structure T>q exactly 
defines the regular reduced Dirac structure on Mq. 

5 Singular dynamics 

In this section a Hamiltonian formalism is described corresponding to the singular reduced Dirac 
structure T>q of Definition 1181 This formalism defines the dynamics corresponding to an implicit 
Hamiltonian system (Mo, T>q, Ho) on the topological space Mq. We show that if (Mo, T>q, Hq) is the 
reduction of the implicit Hamiltonian system (M, T>, H) to Mo, then the G-projectable solutions of 
(M,T>,H) project to solutions of the reduced system (Mq,Vq,Hq). 

First let us define a Hamiltonian formalism on a topological space in the spirit of Sikorski differ- 
ential spaces (see |2E1^11)- Consider a topological space Mq together with a subalgebra C°°(Mq) 
of the continuous functions on Mq, called the set of smooth functions on Mo- A continuous curve 
j(t) on Mo is said to be smooth (see |26j ) if /o o 7 is smooth, as a function from (a subinterval of) 
M. to R, for every Jq G C°°(Mo). Let X denote a derivation on C co (Mq). An integral curve of X 
through some point xq G Mo is defined (see UHl) as a smooth curve 7(t) for which, cf. (|4.6|) . 

^/o(7(*)) = X[M(l(t)), Vt, V/o G C°°(M ), 7 (0) = x . (5.1) 

Let T>q be a topological Dirac structure on Mo, consisting of a generalized Poisson bracket {•, -}o : 
C°°(Mo) x C°°(Mo) — > C°°(Mo) (that is, the Jacobi identity does not necessarily hold) and a vector 
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space A of derivations on C°°(Mo). Furthermore let Hq G C°°(Mq) be a smooth function on Mq, 
called the Hamiltonian function. Notice that {-,Hq}q : C°°(Mq) — > C°°{Mq) defines a derivation on 
C°°(Mq) by {;H } [f } := {fo,H } , f G C°°(M ). Furthermore, if X is a derivation on C°°(M ) 
and x G M , then : C°°(M ) -> M is defined by (X(x))[/ ] := X[/ ](x), /o G C°°(M ). The 

three-tuple (Mq,T>q, Hq) defines an implicit Hamiltonian system in the following way: 

Definition 19. A smooth curve j(t) on Mq is called an integral curve (or, solution,) of 
(Mq,T>q, Hq) if there exists a derivation X on C°°(Mq) such that j(t) is an integral curve of 
X and 

X( 7 (t)) - {, ff }o(7(*)) e A( 7 (i)), Vt, (5.2) 
Z[iT ](7(t)) = °> v ^ v ^ G ^. (5.3) 
The implicit Hamiltonian system (Mq,T>q, Hq) is defined as the total set of integral curves j(t) of 

(Mq,Vq,Hq). 

If Mq is a smooth manifold, Definition El of an implicit Hamiltonian system equals Definition 
El given in Section [2] (with T>q defined by (|4.12|l ). However, since in general Mq is not a smooth 
manifold but only a topological space, the implicit Hamiltonian system (Mq,T>q, Hq) cannot be 
written as a set of differential and algebraic equations. As for implicit Hamiltonian systems defined 
on manifolds, the implicit Hamiltonian system (Mq,T>q, Hq) is energy conserving, cf. (|2.13|) . 

^(7(*)) = X(Ho)m) = {H ,H } (^t)) = 0, Vt. (5.4) 
Remark 20. Equation (|5.2[) implies that 

J t fo(l(t)) = {fo,H }oMt)), Vt 5 V/ G Av , (5.5) 

where A Vo = {/o G C°°(M ) | Z[/ ] = 0, \/Z G A}. However, (53} does not imply (Q. Even in 
the regular case it is not true that Z being a local section of A is equivalent to Z[fo] = 0, V/o G -4x> . 
A counterexample can easily be constructed (by considering a suitable noninvolutive distribution 
A). ♦ 

Remark 21. If A = (i.e., A = 0, that is, no constraints), then A Vq = C°°(Mq) and I|5.21I5.3|1 are 
equivalent to (|5.5|) . In this case, the Hamiltonian dynamics defined by (|5.5|) is exactly the singular 
reduced Hamiltonian dynamics as defined in ^Ql |H1 ESI 123 l2Zj ■ ♦ 

Recall that implicit Hamiltonian systems defined on manifolds define a set of differential and 
algebraic equations and, as a consequence, the standard results on existence and uniqueness of 
solutions for ordinary differential equations do not apply. As explained in Section 12 in general 
one cannot expect neither global existence nor uniqueness of solutions of these systems. Therefore 
one cannot expect global existence and uniqueness of solutions of implicit Hamiltonian systems 
on topological spaces as defined in Definition 1191 In particular, all solutions necessarily lie in the 
constraint space 

M C = {x G Mq I Z[Hq](x) = 0, VZ G A} (5.6) 



19 



(a topological space whose topology is induced from Mo). However, if (Mq,T>q, Hq) is the sin- 
gular reduction of an implicit Hamiltonian system (M,T>, H), then we shall show next that the 
G-projectable solutions of (M, V, H) (if they exist) project to solutions of the reduced system 
(M ,V ,H ). 

Let x(t) be a solution of (M,V,H) with x(0) G N = P' 1 ^). Then by (JUS}, the curve x(t) is 
contained in N. Now assume that x(t) is a G-projectable solution, that is, there exists a projectable 
derivation (i.e., vector field) X on C°°(M) , which projects to a well defined derivation X on 
C°°(Mo), such that x(t) is an integral curve of X (i.e., x(t) = X(x(t))). By (|2~8|) and (f2~l"2j) 

X(z(t)) - {;H}(x(t)) =: Z(x(t)) G A{x(t)), Vt, (5.7) 

F[iT](a^)) = 0, Vt, VY G r ioc (A). (5.8) 

Let Mq be the singular reduced space and T>q be the singular reduced Dirac structure on Mq. Since 
H is assumed to be G-invariant, its restriction to N projects to a well defined function Hq G C°°{Mq) 
defined by Hq on = H\n; for the definition of C°°(Mq) see Definition ^] in Section 0J Define the 
singular reduced implicit Hamiltonian system (Mo, T>q, Hq) as in Definition 1191 Project the curve 
x(t) to Mo to obtain the smooth curve 7(t) = w(x(t)) on Mo- Then ^(t) is an integral curve of 
the derivation X. Indeed, take an arbitrary Jq G C°°(Mo) and let / G C°°(M) G be such that 
/o°t = f\N- Then 

J t foh(t)) = j t f(x(t)) = X[f](x(t)) = X[/o]( 7 (t)), W, (5.9) 

where we used the fact that x(t) is an integral curve of X, cf. I|4.6|) . an d that X projects to a 
derivation Jt on C°°(Mo). Furthermore, if Yi, . . . ,Y m is a basis of projectable local sections of A, 
projecting to a basis Yi, . . . , Y m of A, then it follows that 

= Yj[H](x(t)) = Yj[fZo](7(t)), Vt, j = 1, . . . ,m, (5.10) 

which yields equation (|5.3|) . It remains to be proved that ()5.2j) is satisfied. Notice that by (|4.3j) the 
derivation {-,H} projects to a well defined derivation {-,Hq}q on C°°(Mo). Since also X projects 
to a derivation X on C°°(Mo) it follows that the derivation Z projects to a well defined derivation 
Z on C°°(Mo). Since Y±, . . . , Y m is a projectable basis of local sections of A it follows from ()5.7|) 
that at each point xq on the curve x(t), one has 

Z(xq) = c 1 Y 1 (x ) H h c m y m (x ), (5.11) 

for some constants c\, . . . , c m G R. We claim that 

Z(-y ) = ciFi(7o) + • • • + c m r(7o), 7o = 7r(x ). (5-12) 
Indeed, take an arbitrary f G C°°(M ) and let / G C°°(M) G be such that / o vr = /|at. Then 

(Z(7o))[/ ] = Z[/o](7o) = ^[/](x ) 

= (c 1 F 1 [/] + --- + c m F m [/])(x ) 

= (c 1 y 1 [/o] + --- + c m y m [/ ]) ( 7o ) 

= (ciFi(7o) + • • • + c m y m ( 7 o)) [/o], (5.13) 
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which yields (|5.12|) . Since Yi,...,Y m forms a basis of A it follows that Z{^{t)) G A(7(t)), Vt. 
Therefore (|5.2|) is satisfied, which implies that ~f(t) is an integral curve of the reduced implicit 
Hamiltonian system (Mq,T>q, Hq). We have obtained the following 

Proposition 22. Every G-projectable solution x(t) of (M,D, H) with x(0) G P -1 (0) projects to a 
solution "f(t) = Tr(x(t)) of the singular reduced implicit Hamiltonian system (Mq,T>q,Hq). 

Remark 23. Suppose that the implicit Hamiltonian system (M, T>, H) is of index 1. As remarked in 
Section[21 the system can be restricted to an explicit Hamiltonian system on the constraint manifold 
M c defined by a generalized Poisson bracket denoted by {•, -} c : G°°(M C ) x C°°(M C ) -> G°°(M C ). 
The G-action leaves the manifold M c invariant and it follows that G is a symmetry Lie group of the 
explicit Hamiltonian system on M c , i.e., L^ Mc {f,g} c = {L^ M J,g} c + {f,L^ M j} c , Vf,g G C°°(M C ) 
and L^ M H c = (where H c = H\m c ), for all £ G Q. The corresponding equivariant momentum 
map is given by the restriction of P to M c . Furthermore, every solution of the (restricted) system 
is G-projectable [2S1 03 |H] ■ We can use the singular reduction theory developed in JTU1 IHl 12*31 [2*7] . 
or equivalently the theory developed in this paper by considering A = 0, to reduce the system 
to a Hamiltonian system on the singular reduced space (M c )q. The reduced generalized Poisson 
bracket ({•, -}c)o is defined analogously to l)4.3|) . The reduced dynamics is given by equation 1)5.5)1 . 
with {•, -}o replaced by ({•, -}c)o an d *4/p = C°°((M c )o). Global existence of solutions now follows 
from Proposition 1221 Furthermore, if the G-action is proper then also uniqueness of solutions of 
the singular reduced system can be proved Ell • ♦ 

6 Orbit type decomposition 

Consider a symplectic manifold (M,uj) admitting a symmetry Lie group with a corresponding 
momentum map. Let Mq denote the singular reduced space and {•, -}o the singular reduced Poisson 
bracket. The singular reduced Hamiltonian dynamics is defined by equation 1)5.5)) . cf. B,emark|211 In 
[U 03 El ini 121 E3 it is shown that the space Mq may be decomposed into a family of symplectic 
manifolds, called pieces. The decomposition is by orbit type and defines a stratification of the 
singular reduced space Mq. Furthermore, the Hamiltonian flow corresponding to 1)5.5)1 leaves the 
pieces invariant and restricts to a regular Hamiltonian flow on each of the pieces. In this section 
we show that these results can be generalized to singular reduced implicit Hamiltonian systems. 
We will treat only the special class of Dirac structures for which the generalized Poisson bracket is 
nondegenerate. 

Consider an implicit Hamiltonian system (M, T>, H) with a Dirac structure T> as defined in (|2.8[) . 
or 1)4.1)) . where the generalized Poisson bracket is assumed to be nondegenerate and defined by a 
nondegenerate two-form uj on M, i.e., the generalized Poisson structure J : T*M — > TM is given 
by uj- 1 , that is, {f,h} = u(X f ,X h ), f,h £ G°°(M), where X f is defined by df = u(X fr ) and 
analogously for X^. Notice that we do not assume that a; is a closed two- form. This means that 
{■, •} does not necessarily satisfy the Jacobi identity. 

Let G be a symmetry Lie group of (M,V,H) as in Proposition From now on we will assume 
that the action of G is proper. Assume that there exists an AcT-equivariant momentum map 
P : M -> q* such that 

dPj: = u(£ M , ■) (equivalently £ M = JdP ( ) and dP c G r(A°), V£ G g. (6.1) 
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Notice that this implies, but is not equivalent to, (|3.7|) . 

Remark 24. Symmetry groups as described above commonly occur within the class of constrained 
mechanical systems as described in Example Consider the distribution A = ker^4 T (g) on Q and 
let G be a Lie group acting properly on the configuration space Q leaving A invariant, that is, 
L^ Q Ti oc (A) C r ioc (A), V£ G g. The action of G on Q lifts to an action of G on M = T*Q as follows 
(recall that u = dp A dq denotes the canonical symplectic form on M): define the infinitesimal 
generator £m to be the Hamiltonian vector field corresponding to the function P^(q,p) = p T £,Q(q), 
i.e., dP^ = lo(^m, ^£ ^ 0- Now recall that A is defined as the image of the matrix (|2.19fl . Then 
by construction, L^ h[ uj = and L^ M Ti oc (A) C r^ oc (A), V£ G g, and therefore G defines a symmetry 
group of the implicit Hamiltonian system on M. The ^4d*-equivariant momentum map is defined 
by P(q,p)(£,) = p T £,Q(q), V£ G g, and hence satisfies the first condition in (|6.1|) . If we assume 
furthermore that the symmetry group is "horizontal", i.e., £q G keiA T (q), V£ G g, then also the 
second condition in (|6.1|) is satisfied. ♦ 

The manifold M can be decomposed into submanifolds as follows jHlEUE^j. Let K be a compact 
subgroup of G and define Mrx) to be the set of points in M whose stabilizer group G x = {g € G \ 
4>(g,x) = x} is conjugate to K, i.e., 

M {K) = {x G M \3g £G such that gG x g^ = K}. (6.2) 

Notice that since the G-action is assumed to be proper, every stabilizer group G x , x G M, is a 
compact subgroup of G. M^) is a submanifold of M called the manifold of orbit type (K). On the 
set of compact subgroups of G define an equivalence relation by saying that K ~ K if and only if 
K is conjugate to K. The equivalence class of K is denoted by (K). As (K) runs over the set of 
equivalence classes, the manifolds M^ K ^ partition M. Since the G-action is proper this partition is 
locally finite. This is called the orbit type decomposition of M. 

Next we show that the image of the tangent of the momentum map at the point x G M is equal 
to the annihilator in g* of the Lie algebra of the stabilizer group G x , i.e., 

lmT x P = g°, G M, (6.3) 

see also EEU EEH1 EU • Indeed 

?6S^ £m(x) = o cLP e (x) = & (T x P(v))£ = 0, Vi> G T X M f G (Im T a .P)°, (6.4) 

where we used (j6.1|) and the fact that u is nondegenerate. This yields that g^ = (Im T X P)°. Taking 
the annihilator of both sides yields (|6,3|) . 

Equation (|6.3|) implies that the tangent of the restriction of the momentum map P to the manifold 
M(jn has constant rank equal to the codimension of K in G. It follows that the intersection 
P _1 (0) n Mt^\ is a smooth submanifold of M. Furthermore, the manifold P _1 (0) D Mr K \ is 
invariant under the action of G. It turns out that the quotient (Mo)fK) := (P 1 ^) H M^)/G = 
7r(P _1 (0) D Mix)) is a smooth manifold [2l|2lll23|. Consequently, the singular reduced space Mq 
is decomposed into a disjoint set of manifolds, called pieces, 

Mo = I] (Mo) w , (6.5) 
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where (K) runs over the set of conjugacy classes of compact subgroups of G. Since the orbit type 
decomposition of M is locally finite, the decomposition of Mq is also locally finite. 

Next let us define a generalized Poisson bracket on each of the manifolds (Mo)nn. For clarity of 
exposition, consider the following commuting diagram: 

M 




(o) n M (K) 



7T 



7T 



(K) 



M 



(M )(K) 



Here irj^) denotes the inclusion map and t^ik) the restriction of tt to P 1 (0) DMr^)- The inclusions 
and i® K j are defined such that the diagram commutes. Define the set of Whitney smooth 



Hk) 



functions W°° ((Mq) ^) on (Mo) on as follows: a continuous function /q on (Mo) on is said to 
be a Whitney smooth function if there exists a smooth G-invariant function / G G°°(M) such 
that /o o 7T(^) = /|p-i(o)nAf (iC ) ■ I n f ac t) W°°((Md)(.k")) is equal to the set of functions obtained by 
the restriction of the functions in G°°(Mo) to (Mo) on (this is why it is called the set of Whitney 
smooth functions). Indeed, the G-invariant function / descends to a smooth function /o on Mq, 
whose restriction to [Mq)i k \ is precisely fo. This can be seen as follows: 



fo ° i 



fo 



O TT O L 



(K) 



f ° t 



o / 



(JO 



/ 



o / 



fo 



O 7T 



(if), 



(6.6) 



and since 7ron 1S surjective the result follows. Analogously to (|4.3j) define a generalized Poisson 
bracket on W°°((M ) iK) ) by 



{fo,ho}( K ) 07T (K) = {/, ft}[p-i(o)nM ( 



A') 



(6.7) 



We need to show that the bracket is well defined and does not depend on the choice of G-invariant 
extensions / and h. If we can prove that for every / G G°°(M) G the flow of the Hamiltonian 
vector field X f = {•,/} preserves the submanifold P _1 (0) n M (K) , then 1 = {/ G G°°(M) G | 
7|p- 1 (o)nM (Jf) = 0} is a Poisson ideal of C°°(M) G and it follows that the bracket is well defined, 
see e.g. 0, Appendix B, (5.2) & (5.3). Now consider / G G 00 (M) G , that is, </>*f = f, for all g G G, 
where (f) g : M —> M denotes the diffeomorphism given by the action of the element g G G on M. 
Thus, since the action is symplectic, it follows that 4>* g Xf = X^f = Xf for all g G G, which implies 
that the flow tp{ of commutes with <fi g for every g G G. We need to show that V/ preserves 
M( K y Since Mon = G ■ Mr, where Mk = {x G M \ G x = K} (which is a submanifold of M) and 
the flow of Xf and the G-action commute, it is enough to show that ip{ preserves Mk- Now, for 
every g G K and x G Mx 

<P g ^{(x)) = = i/>!(x), (6.8) 

since 4> g {x) = x. Therefore, K C G^f^ (the stabilizer group of tj){(x)). Suppose g G G^j^y i.e., 
4> g (ip{ \x)) = tp{(x). Since the flow and the G-action commute, this implies that ijj{(4> g (x)) = ip((x). 
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Applying \jj_ t to this relation yields (f> g (x) = x so g G K, i.e., G^f^ C K. It follows that 
^tf> f (x) = ^ anc ^ therefore ip{ (x) G M/<-. So we have shown that for any G- invariant function / 
the flow preserves the submanifold P _1 (0) D Mix)- This implies that the bracket H6.7|) is well 
defined. 

Remark 25. The set of Whitney smooth functions W°°((Mo)(k)) is dense in the set C°°((MoW)) 
of smooth functions as defined by the differential structure on (MoV#-)- Indeed, the pullback to 
P _1 (0) fl M(x) of a smooth function /q G C 00 '((Mo) ik)) compactly supported on (Mo)(^) can be 
extended to a smooth G-invariant function / on M. Thus the bracket in 1)6 .7)1 is a well defined 
generalized Poisson bracket {•, -}(k) on C°°((Mo)(in). The associated generalized Poisson structure 
J(K) is the inverse of the nondegenerate two-form uJ(k) defined by the condition that the pullback 
of uj {k) to P -1 (0) n M(jf) equals the restriction of uj to P -1 (0) nM^, see also E31 E3j - ♦ 

Having defined a generalized Poisson bracket on (Mo)on we now show that the distribution A 
projects to a distribution (A)on on the manifold (Mq)(x)- Recall from Section |I] that there exists 
a basis Xi, . . . ,X m of local sections of A such that [£jw,Xj] G 5J/ O o V£ G g, j = 1, . . . , m. By 
Proposition El it follows that 

[Cm, X,] G T ioc (A nn V£ G fl, j = 1, . . . , m. (6.9) 

If we make the assumption that A n V = 0, then it follows that [£mi XA = 0, V£ G g, j = 1, . . . , m. 
This means that the flow of Xj commutes with the G-action and therefore preserves the submanifold 
P _1 (0)nM(^). This implies that Xj is tangent to P'^OjnM^j and restricts to a local vector field 
Xj on P _1 (0) fl Mrx)i f° r an 3 = 1) • • • j m - The basis Xi, . . . , l m projects to a set of independent 
local vector fields Xi, . . . ,X m on (Mq)(x) which define a distribution (AV#n on (Mq)(x)- 

Remark 26. If the system is of index 1, then it is always true that A(x) n V(x) = 0, \/x G M c 
(see [Hlini)- ♦ 

The generalized Poisson bracket {•,•}(_&-), with its generalized Poisson structure J (%;)■> an d the 
distribution (A)^ define a Dirac structure *D{k) on (Mo)(#), given in terms of its local sections 
by 

{(X,a) G X loc {(M ) (K) ) x nj oc ((M ) w ) | 

x - j (K) a g r ioc ((A) (Jf) ) , a g r ioc (((A) (K) )°)} . (6.10) 

The results in Section pimply that X>(jn is exactly the regular reduced Dirac structure on (Mo)(jn. 
We summarize the previous discussion in the following statement. 

Proposition 27. Consider an implicit Hamiltonian system (M, T>, H) with a Dirac structure T> of 

the type given in \2. 8\) . defined by a nondegenerate generalized Poisson bracket {-, •} and a smooth 

vector subbundle A of TM . Assume the system admits a symmetry Lie group G acting properly 

on M and satisfying the hypotheses of Proposition^ Suppose also that the action admits an Ad*- 

equivariant momentum map P satisfying hb\l\) . Assume furthermore that A n V = and that 

the assumptions in Propositions^^ and \l$ are satisfied. 8 Then the singular reduced space Mo is 

8 Notice that A and V being constant dimensional distributions and A n V = imply that A + V is a constant 
dimensional distribution, i.e., a smooth vector subbundle of TM. 
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decomposed into a disjoint set of manifolds (Mo)(in, cf. h6.5\) . called pieces. The system (M,T>,H) 
reduces, by means of regular reduction, to an implicit Hamiltonian system ( (Mo )(k)i ^(k)i H{k)) 
on each of the pieces. Here T>r^\ is the Dime structure defined in \6.1b)) and Hr^\ G C°°((Mo)(x)) 
is defined by H {K) o tt {k) = H\p-i {0)nM{K) . 

Remark 28. The only reason why we assumed that the Dirac structure T> is defined by a non- 
degenerate generalized Poisson bracket and the momentum map P satisfies (|6.1|) . is to show the 
equality (|6,3[) . This equality implies that the subsets P -1 (0) fl M(jn are smooth manifolds and 
consequently the pieces (Mq)[k) are smooth manifolds as well. The construction of the reduced 
Dirac structure T^(k) however is completely general. This means that the result of Proposition is 
valid for general Dirac structures and Ad*-equivariant momentum maps of the form (|3.7|) . as long 
as it is certain that the pieces (Mq)(x) are smooth manifolds. The same will hold for the results in 
the remaining part of this section. 9 ♦ 

Next we show that the regular reduced implicit Hamiltonian systems ((Mo)(x),XW), -H(in) are 
precisely the restriction of the singular reduced implicit Hamiltonian system {Mq,T>q,Hq) to the 
pieces (Mo)^ K y Indeed, first we show that the inclusion : (Mq)^ — > Mq is a Poisson map. 
Consider f ,h G C°°(M ) and their restrictions f Q ,h G W°°{{Mq) {k) ) and let f,h G C°°(M) G 
be such that fo o ir = /|p-i(q), and analogously for h. The definition of the brackets and the 
commutativity of the diagram implies that (notice that by (Jfi.fijl /o o h^k) = /|p- 1 (o)nAf (Jf) and 
analogously for ho) 

{fo,h }(K)°^(K) = {f,h}oL {K) = {f,h}oLol {K) = {/ ,^o}o°vroi (x) = {/ , /io}o°4W) 07r (K)- (6-11) 

Since tt(k) is surjective it follows that 

{fo,ho}(K) = {/o>o}o|(A/ ) w - (6.12) 

By construction, it is immediately clear that (A)nn is the restriction of A to W°° ((Mq) r^) , defining 
by denseness a distribution on C°°((Mq)ij^\). We conclude that the regular reduced Dirac structure 
D(^) is precisely the restriction of the singular reduced Dirac structure Vq to the piece (M )(x)- 
Next, consider a solution 7 (i) of (Mq,T>q, Hq), as defined in Definition fTTH Since X\ , . . . , X m is a 
basis of local sections of A there exist local functions Cg, . . . , c™ G C°°(M ) such that 

X(j(t)) - flb}o(7(*)) = cS( 7 (t))*i(7(*)) + • • • + c rt ( 7 (t))l m (7(i)) G A( 7 (t)), Vt. (6.13) 
Consider the derivation 

y = {•, # }o + clXt + ... + c^X m (6.14) 

on C°°(M ). This is the projection of a local vector field Y = {■, H} + c 1 X 1 + ... + c m X m on M, 
where c J G C 00 (M) G are such that Cq o 7r = c' lp-im) and {Xi, . . . ,X m } denotes the projectable 
local basis of A. Since the flow of this vector field commutes with the G-action, it preserves the 
submanifold P~ 1 (0) fl M^) (following the argument above Remark I25|). It follows that the flow 
corresponding to the integral curve j(t) preserves the pieces (Mq)(k) and therefore 7 (i) restricts to 

9 One could even think of further generalizing the results by allowing the pieces to be topological spaces instead 
of manifolds. 
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a smooth curve 7(t) on (Mo)rjn. The vector field Y restricts to a vector field Y' on (Mo)(x)- By 
construction it follows that 

Y' = {-, H (K) } {K) + c l X[ + ... + c^X' m , (6.15) 

where Cq = Cg|(AfoW) an d Xj is the restriction of Xj to (M )(x)> j = 1, . . . ,m. The curve 7 (i) is 
an integral curve of Y' . Indeed, 

P° m) = Jt Ul{t)) = 1[/o](7(t)) = ^[/o](-y(*)) = ^'[/o](7(*)) ; V/o G W°°((M ) W ), (6.16) 
and since W°°((Mq)(x\) is dense in C°° ((Mq) (jq) the result follows. Furthermore, it is clear that 

^(7(*))- j( K )m))dH {K) m)) = Y'^(t))-{;H {K) } {K) m)) G r, oc ((A) w ) ( 7 w), vt, (6.17) 

and 

Z'[F w ]( 7 (i)) = 0, Vt, VZ' € r ioc ((A) w ) , (6.18) 

which means that 7(i) is a solution of the regular reduced Hamiltonian system ((Mo) on , 'D{k\ , H^). 
In conclusion we have proved the following: 

Proposition 29. Consider the conditions in Proposition \21\ The regular reduced implicit Hamil- 
tonian systems ( (Mp )(#•), Z> on > -Hon) a? "e exactly the restrictions of the singular reduced implicit 
Hamiltonian system (Mq,T>q, Hq) to the pieces (Mo) on. A solution j(t), with 7(0) € (Mo) on, of 
(Mo, T>q, Hq) preserves the piece (Mo) on and restricts to a solution j(t) of ((Mo)on,^on,-ffon)- 

Finally, remark that since the pieces (Mo) on are smooth manifolds, each implicit Hamiltonian 
system ((Mo) on, T^(K)i H(K)) can be written as a set of differential and algebraic equations (DAE). 
The singular reduced implicit Hamiltonian system (Mq,T>q, Hq) can thus be written as a collection 
of DAEs, one (set of differential and algebraic equations) on each piece. 

7 Conclusions 

In this paper we studied the singular reduction of implicit Hamiltonian systems admitting a sym- 
metry Lie group with a corresponding equivariant momentum map. The results extend the singular 
reduction theory known for explicit symplectic or Poisson Hamiltonian systems [3J 1101 1111 1231 
1241 127j . The main result is a purely topological description of the reduced implicit Hamiltonian 
system using the definition of a topological Dirac structure. In particular, the reduced space is not 
assumed to be a smooth manifold. The dynamics corresponding to this system are defined and it 
is shown that the projectable solutions of the unreduced system project to solutions of the reduced 
system. If the symmetry Lie group acts freely and properly and the value of the momentum map is 
regular, then the singular reduced implicit Hamiltonian system equals the regular reduced implicit 
Hamiltonian system as described in [H] - Finally, under certain conditions, the singular reduced 
space can be decomposed into a set of smooth manifolds called pieces. It is shown that the singular 
reduced implicit Hamiltonian system restricts to regular reduced implicit Hamiltonian systems on 
all the pieces. 
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